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We study the growth of a directed network, in which the growth is constrained by the cost of adding 
hnks to the existing nodes. We propose a new preferential-attachment scheme, in which a new node 
attaches to an existing node i with probability n(fci) oc , where ki is the number of outgoing links 
at i. We calculate the degree distribution for the outgoing links in the asymptotic regime {t oo), 
rife, both analytically and by Monte Carlo simulations. The distribution decays like kfi'^ /r{k) for 
large k, where /i is a constant. We investigate the effect of this preferential-attachment scheme, 
by comparing the results to an equivalent growth model with a degree-independent probability of 
attachment, which gives an exponential outdegree distribution. Also, we relate this mechanism to 
simple food-web models by implementing it in the cascade model. We show that the low-degree 
preferential-attachment mechanism breaks the symmetry between in- and outdegree distributions 
in the cascade model. It also causes a faster decay in the tails of the outdegree distributions for 
both our network growth model and the cascade model. 

PACS numbers: 89.75.Fb, 89. 75. He, 87.23.Cc, 02.10.Ox 



I. INTRODUCTION 



Many real-world networks that have been studied recently, such as the internet, the World Wide Web, and social 
relations have a scale-free structure, whose degree distribution, P{1), follows a power law for large /, the number of 
links connecting to a single node. A scale-free network can be constructed using a preferential-attachment mechanism 
in which a new node attaches to an existing node, i, which has li links with a probability = h/^j Ij !]■ This 
form of preferential attachment means that the highly connected nodes attract more new nodes than the others, 
and it explains the structure of networks like the internet, in which a connection is favorable for both connected 
nodes, and where links arc undirected. Growth according to this mechanism leads to self-organization into a scale-free 
structure X 2> 1, i., ^] ■ 

In this paper, we propose a new preferential-attachment scheme for growing directed networks |a, in which new 
nodes prefer to attach to existing nodes with lower degree. Such a mechanism could play a role in some transportation 
networks |i8„^] such as food webs (10., Jd,^, 13] or power-grids tl^il^], which can be considered directed. The hnks 
in such networks transport some sort of resource to "feed" the nodes. Each node gets resources only by feeding on 
another one. (The initial nodes can be considered as sources.) We take the direction of the resource flow as the 
direction of the link. Therefore, an outgoing link transports resources from a "prey" node to a "predator" node. If 
there is more than one predator that feeds on a prey node, then the resources supplied by the prey have to be shared 
by the predators. This implies that for a conventional growth model, the attractiveness of a node in such a network 
should increase with the number of its incoming links (indegree) , and decrease with the number of its outgoing links 
(outdegree). In other words, the cost of adding links to an existing node should increase due to limited capacity, as 
the number of nodes that feed on it increases. This cost puts a constraint on the growth of the network A good 
example of such a system, the world-wide airport network, is studied in Ref. [T^ . 

Therefore, we propose a new growth scheme, in which a new node attaches to an existing node, i, with probability 
Il{k'^,ki) (X (k'^/ki)'^, where ki and k[ are the number of outgoing and incoming links at node i, respectively. This 
means that new nodes prefer connecting to existing nodes with more incoming and fewer outgoing links because they 
provide more resources per outgoing link. Here we study a special case in which k[ = m for all nodes. 

As mentioned above, we believe that such a mechanism could play a role in real systems such as food webs, which 
provide a good basis for testing our predictions. There has been significant progress in this field in the last decade. 
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FIG. 1: Illustration of the growth with A'o = 4 and m — 3. Nodes 1, 2, 3, and 4 are the initial (source) nodes connected to 
unspecified nodes outside the network. Nodes 5, 6, and 7 are the nodes added at timesteps 1, 2, and 3, respectively. See text 
for details. 

especially in development and understanding of static models, which have been quite successful at reproducing webs 
with structural properties in agreement with the empirical data [ToLllTlll3 .ri3. 17. 18. 19, 20] . Despite their success, 
most, if not all, proposed static food-web models in the literature such as the cascade model |2l[ and the niche 
model [T^ are phenomenological models. They do not employ physical principles in the assembly process — or at 
least, the rules they employ do not yet have any physical interpretations. Although the mechanism we propose in 
this paper is quite naive, we believe that a static food-web model should employ such rules, which can be related to 
physical principles. 

The structure of the paper is as follows. We describe the model, derive the outdegree distribution, and compare the 
results with Monte Carlo simulations in Sec^l Some details of the calculations are given in Appendix El In Sec lIIII 
we compare the proposed model to an equivalent model with a degree-independent probability of attachment. The 
derivation of the outdegree distribution for this model is provided in Appendix [CI In Appendices ^ and (CJ we show 
that the outdegree distributions for these growth models can be scaled onto a single scaling function for large m. In 
Sec mil we also relate our results to some food- web models. We summarize our results in Sec. lIVI 



The growth begins with A'o > m nodes with m incoming links, each connected to unspecified nodes outside the 
network. These initial nodes can be considered as sources, and the only purpose of their links is to make sure that 
each node has the same number of incoming links. We add one node with m incoming links at each time step, which 
connects to m different nodes. The direction of the new links are the same as the direction of the resource flow: from 
the existing node to the new one. An example of such a network after the first few time steps is shown in Fig. 

As we stated above, an existing node acquires a new link with probability Il{k[,ki) cx (k'Jki)'^. However, for the 
case in which all nodes have m incoming links, /c^ = m is just a scale factor. Without changing the proportionality 
for fci 1, we replace ki by fc^ + 1 in order to keep Il{ki) finite for ki = 0. Here, we study a special case of this 
attachment probability, H(fci) oc k^'^ , in which 7=1. Therefore, the probability of attachment is given by 



II. MODEL AND RESULTS 
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Here, Nk is the number of nodes with outdegree fc, and Nt is the total number of nodes. We note that Nk, and 
therefore, Zm, depends on time. 

We use the rate-equation approximation of Refs. ^,|^. In the hmit of Nt ^ 1, the rate equations for the outgoing 
hnks are: 

^ - 1 - mn(0)7Vo , (3) 

= mIV{k - l)Nk-i - mn{k)Nk for fc > . (4) 

Substituting nkNq^ for N^, and using dN^/dt — 1, yields corresponding equations for the density of nodes: 

dno I — uq 



dt Nt 



mU{0)no , (5) 



^ = mU{k - l)nk^i ^ Uk [N^^ + mU{k)] for fc > . (6) 

We are interested only in the asymptotic regime (t — > oo), in which duk/dt = 0, and the distribution reaches a steady 
state. The steady-state solutions of Eqs. (jSJ and © are 

^ (7) 



k[{k + + z*^/m 
where 



1 + m/z. 



'''^ for fc> , (8) 



z;, = lim Z„jNrr = V ^ . (9) 

1=0 



The solution of the recursion relation Eq. (|SJ) can be obtained by a few straightforward algebraic manipulations. First, 
we convert the right-hand side of Eq. (jSJ into a product, using Eq. ((T)): 



im/z*J'^{k + l) 



k 

n I * for fc > . (10) 

J-A 7 + l + m z* 



1 + m/z*„ -11 J + 1 + m/ 

Then, using the identity r(a; -I- 1) = xT{x)^ we write the product in terms of a ratio of Gamma functions to obtain 
the final form of the outdegree distribution: 

nl = [k + l){m/z:,f ^ I ' y' forfc>0. (11) 

We note that Eq. I|ll() holds for fc = as well. (The outdegree distribution, Eq. (|ll|l . is analogous to the predator 
distribution in a food web. We use the terms in- and outdegree distribution interchangeably with prey and predator 
distribution, respectively.) The outdegree distribution, n^, is shown in Fig. [3 for m — 1,3, and 5. The method we 
use to obtain zj^ is explained in Appendix 1X1 Also, in Appendix ^ we prove that the scaled distributions mn^ vs 
u = k/m collapse onto a single scaling function, 

p{u) = ^eM~^uyA), (12) 

for large values of m. 

We also check our theoretical results by comparing them to results of Monte Carlo simulations. The growth scheme 
in the simulations is implemented as follows. We begin the simulation with Nq > m nodes with m incoming links 
each. At each time step, we add one new node. We pick a random existing node from the network, calculate the 
probability of attachment to that node, and generate a random number. If the random number is smaller than the 
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FIG. 2: (Color online) Outdegree (predator) distributions, nl, for the growth model with low-degree preferential attachment, 
for m = 1, 3, and 5 with A'o ~ 10. The results of the Monte Carlo simulations (the symbols *, x, and +) are in excellent 
agreement with the theory (the symbols o, □, and o). As fc is a discrete variable, the lines connecting the symbols are merely 
guides to the eye. The simulations were performed up to a system size of A'r = A^o + 10^ nodes, and averaged over eight 
independent runs. Inset: The same distributions shown on a log-linear scale. The tails decay faster than an exponential, 
indicating that the topology of the network is relatively homogeneous. 
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FIG. 3: (Color online) Outdegree (predator) distributions, nl, for the equivalent growth model with degree-independent 
probability of attachment, for m = 1, 3, and 5 with A'o = 10. The results of the Monte Carlo simulations (the symbols *, x, 
and +) are in excellent agreement with the theory, which predicts exponential decay (the symbols o, □, and o). Compare to the 
degree-dependent case in Fig. |21 The simulations were performed up to a system size of iVT = iVo + 10*^ nodes, and averaged 
over eight independent runs. Inset: The same distributions shown on a log-linear scale along with the theoretical results for 
the degree-dependent case (gray curves - turquoise online) in Fig. |21for comparison. 



attachment probability, we create a link between the new node and the randomly picked node. If not, we randomly 
pick another node from the network and repeat the process until the new node has m links. (We do not pick the 
same node twice.) The direction of a new link is from the existing node to the new one. The simulation stops when 
the system size reaches iVmax + -^o nodes. We take A'^ax — 10^ and average the results over eight independent runs. 

The results of the Monte Carlo simulations are also shown in Fig. [21 They are in excellent agreement with the 
theory. The deviations in the last three nonzero data points are consistent with the number of runs used for averaging. 

In order to clarify the role of the low-degree preferential attachment, we provide some comparisons with similar 
models in the next section. 
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FIG. 4: (Color online) Degree distributions for the cascade model with and without low-degree preferential attachment with 
d= 10 and S = lO". The error bars for the Monte Carlo data (averaged over eight runs) are smaller than the symbols. "Total" 
represents the distribution for the total number of links, k + k' . (a) The distributions for the original cascade model, in which 
the probability of attachment is x — d/S (the symbols, o, □ and v)- The in- and outdegree distributions obtained from the 
Monte Carlo simulations practically overlap in agreement with the theory, Eq. 11411 (black line), (b) The distributions for the 
cascade model with low-degree preferential attachment with x{ki) = d/\S{ki + 1)]. The in-outdegree symmetry of the original 
cascade model is broken in the new distributions. See text for details. 
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FIG. 5: (Color online) Scaled in- and outdegree distributions for the cascade model with (the symbols, A and o) and without 
(solid circles) low-degree preferential attachment. In the axis titles, n„ represents the degree distributions, n^, n^/, and 
p{in=oiit)^^-j^ and (k) represents their respective mean values, which are 5 for the original cascade model, and 2.36 for the 
modified model. Each distribution is scaled using its respective (k). The in- and outdegree distributions have the same mean 
for each case. Inset: The same distributions shown on a log-linear scale for comparison. 



III. THE ROLE OF THE LOW-DEGREE PREFERENTIAL ATTACHMENT AND RELEVANCE TO 

SIMPLE FOOD- WEB MODELS 

A. A growing network with degree-independent probability of attachment 

We first compare our model to an equivalent model with a degree- independent probability of attachment, IT = 1/Nt- 
In this case, 11 is the same for all nodes, however, it is time dependent. As we prove in Appendix El a degree- 
independent n yields a network with an exponential outdegree distribution as seen in Fig. |21 This shows that the 
low-degree preferential attachment is responsible for the peaked form of n^, and also for the faster decay in the tail. 
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B. The cascade model 



Secondly, we test the effect of the low-degree preferential attachment in the cascade model |20|, |21| , which is a simple 
wiring-based food-web model. Although the cascade model has now mostly been replaced by improved variants like 
the niche model it serves our purpose better, because it is quite simple and transparent. The cascade model 
builds a web using a few simple rules as follows. The species are numbered from 1 to S, where S is the total number 
of species. Each species, i, is allowed to prey on a species j, where i > j, with a constant probability x = d/ S, where 
d < S' is a constant. The species with no predators constitute the top species, and the species with no prey constitute 
the basal species The average number of prey (or predator) per species is equal to {d/2){l — l/S*) k, d/2 for 

large S. Therefore, d/2 is analogous to m in our growth model. 

The calculation of the outdegree (predator) distribution of the webs produced by the cascade model is quite 
straightforward. First, we find the probability of finding a node, i, with k outgoing links: 

V'f'''\k)^(^^'^{\~xf''-''x'' wi\hk<S-i. (13) 
Then we obtain the outdegree distribution by summing over i, and normalizing by 5* |24j: 



p(out) (fc) = i ^ T'f™*^ {k) with k<S . (14) 



The cascade model produces webs with identical in- and outdegree distributions. One can prove this by substituting 
Vt%,{k) for Vt'ik) in p('")(fc) = Etfe+i vt\k)IS. 

In order to see the effect of the low-degree preferential attachment, we modify the cascade model by changing the 
probability of attachment from a constant, x — d/S, to a decreasing function of the outdegree of the prey node, 
x{ki) — d/[S{ki + 1)]. (We note that for both original and modified cascade models, the probability of attachment 
does not depend on time, unlike in our proposed growth model.) Obviously, this modification does not conserve the 
expected total number of links. Nevertheless, it still demonstrates the effect of the low-degree preferential attachment 
in the cascade model. We also point out that the degree-dependent nature of this mechanism could cause different 
degree distributions in the cascade model, depending on the order in the wiring procedure. (See endnote psf.) We 
leave extensive analysis to a further publication and give only these preliminary results here. 

Fig.0]shows the in- and outdegree distributions of the original and modified cascade models for d = 10 and S = 10^. 
One can clearly discern two changes: the new in- and outdegree distributions are not identical, and they are also 
quite different than those of the original cascade model. The new outdegree distribution takes a similar form to the 
outdegree distribution in our model, n|;, in Eq. Hll|l . and the new indegree distribution decays monotonically. For a 
better comparison, we also provide the scaled forms of these distributions in Fig. [S] The figure shows that the tail 
of the outdegree distribution is depressed, while the indegree distribution decays much slower than in the original 
cascade model. All the distributions decay faster than an exponential. The analytical forms of the distributions are 
left for further investigation. 

One could also compare our results with the latest findings in food-web theory. According to Camacho et al. 
[11 El 113, the predator and prey distributions for most empirical food webs obey universal functional forms. The 
universal form of the predator (outdegree) distribution is given by Ppred(fc) = 7(fc + ^,z)/z, where 7(fc,z) is the 
incomplete Gamma function and z is the average connectivity, 2L/S, with L the total number of links in the web. 
This is a continuum approximation to Eq. (|14|l for 5 ^ 1 and a; <C 1. This form of the predator distribution suggested 
by Camacho et al. and the outdegree distribution of our model with low-degree preferential attachment, n^, have some 
similarities, like the cut-off values and Gaussian-like, fast-decaying tails. A detailed assessment, however, requires 
further analysis, which we plan to publish elsewhere. 



IV. SUMMARY 



We have studied an attachment scheme for growing directed networks, in which a new node attaches to an existing 
node i with a probability Tl{ki) oc 1/fci, where ki is the number of outgoing links at i. The motivation behind this idea 
is the following. If the links supply resources to the nodes in the network, and if the nodes can get the resources only 
through other nodes, then a new node will prefer to attach to existing nodes with fewer outgoing links to minimize its 
competitors for resources. We obtained the steady-state outdegree distribution, n^, which decays like kfjl^ /V{k) for 
large k, where is a constant. In order to further understand the effects of the low-degree preferential attachment, 
we also studied an equivalent growth model with a constant probability of attachment, which yields an exponential 
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FIG. 6: (Color online) The m dependence of and the convergence of the ratio of the standard deviation of nl to its mean, 
o/(k), vs m. The graph shows that the standard deviation is proportional to the mean for large m. The ratio a / {k) converges 
to A ~ 0.522 like 1/m. The vs m plot confirms that z^ oc 1/m for large m. 



outdegree distribution. We also implemented this mechanism in the cascade model. The low-degree preferential 
attachment caused a difference between in- and outdegree distributions in the modified cascade model (identical in 
the original model) , as well as a significant depression in the tail of the outdegree distribution and a broadening of the 
tail of the indegree distribution. These results indicate that the main effect of the low-degree preferential attachment 
mechanism is the fast decay in the tails of the outdegree distributions. 

Although we propose this scheme for growing transportation networks, we do not claim that this is the only factor 
that determines the topology of such networks. In fact, we do not believe that growing transportation networks can 
be realistically modeled using a semi-static growth scheme, in which only the new nodes are allowed to make new 
connections while the rest of the network is frozen. Our model also lacks some essential features like the conservation 
of energy or mass. As a result, the networks can grow indefinitely. Nevertheless, the scheme we propose could be 
taken as a first approximation to modeling a growing directed transportation network, in which the probability of 
attachment is a decreasing function of the number of outgoing links. 
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APPENDIX A: NUMERICAL CALCULATION OF £ 



We use numerical methods to determine the values of zj^ for different m. Substituting Eq. into Eq. yields 
a self-consistent equation for z^: 

, ^ {m/z*^fT{\ + m/z;^) 



k=0 



r(fc + 2 + m/z^) 



7 



Asymptotic, Degree-Dep. 
— Asymptotic, Degree-Indep. 

■ — ■ m=10, Degree-Dep. 

■ — ■ m=100, Degree-Dep. 

m=10, Degree-Indep. 

m=100, Degree-Indep. 




FIG. 7: (Color online) Scaled degree distributions (thin curves), mnj! vs fc/m, for the degree-dependent and degree-independent 
models (each for m =10 and 100), shown along with the asymptotic scaling functions (heavy curves), p{u) = ^ exp(— 7rii'^/4) 
and q{u) = exp{—u) with u — k/m. For large values of m, the scaled distributions collapse onto their respective asymptotic 
scaling functions. Scaled distributions are obtained using analytical forms, Eq. IIH and Eq. 1C4II . 



which can also be written as 



<^{l,2 + m/z^;m/z:^) 
1 -I- m/z;; 



(A2) 



where 



where $(1,2 -I- m/ z*^,m/z^) is the degenerate hypergeometric function defined as $(a, 5;z) = J2^=o iWf 
{a)n = a{a + l)(a -I- 2) ... (a -I- n — 1), with (a)o = 1 H^- The right-hand side of Eq. ljA2|l is a smooth function of 
We obtain z^ by solving Eq. 1A2|I numerically for m/z^ For m = 1, 3, and 5, we obtained = 0.6419, 0.3582 
and 0.2467, respectively, which we used to plot the outdegree distribution, Eq. Hll|l . in Fig. [3 The m-dependence of 
zj^ is shown in Fig. 1^1 As seen in the figure, z^ is proportional to m^^ for large m. 



APPENDIX B: THEORETICAL CALCULATION OF THE SCALING FUNCTION FOR THE GROWTH 
MODEL WITH LOW-DEGREE PREFERENTIAL ATTACHMENT 



We scale nl by substituting u — k/m and p{u\m) — mn„iu in Eq- ©■ 



/ I N {mu+ 1) _i 

P[u\m) = p[u-m \m) . (Bl) 

zj:^u[m/z^ + mu + IJ 

By taking the logarithm of both sides, we obtain 

hip[u\m) -\\ip{u-m-^\m) =\n ■ (B2) 

We convert the left-hand side to a derivative, and expand the right-hand side into a series assuming that l/mu, z^u, 
and z^/m — > as m — s- oo: 

19 I z* 

— ^ lnp(w|m) = — - z*-^u . (B3) 

TO ou mu TO 

This equation has a solution of the form 

p{u\m) — Cuexp(— z*Jm + TO7i^/2]) . (B4) 
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We require that p{u\m) should be normaHzable with (u) = 1 (recah that (fc) = m). In the hmit m oo, this requires 
that = a/m, where a is a constant. Using these conditions, we normahze p{u\oo) = p{u) — Cm exp(— 
obtaining C — tt/2 and = 7r/2m. The asymptotic function, p{u), is given in Eq. (|12|l and shown in Fig. [3 The 
excehent agreement with mn'^ vs k/m confirms our resuhs. 

This result also implies that the ratio of the standard deviation of ri^ to its mean, a/ (k), converges to a constant as 
m increases. Therefore, a becomes proportional to m, since (fc) is just m [23|. Fig. El shows that cr/(fc) = A + 0(l/m) 
for large m. We obtain A = limm^oo <^/{k) using p{u): 

A = ^^"'^ 7 ^""^^ = (4/^ - 1)1/2 ^ o_522 . (B5) 



APPENDIX C: THEORETICAL CALCULATION OF OUTDEGREE DISTRIBUTION FOR THE 
GROWTH MODEL WITH DEGREE-INDEPENDENT PROBABILITY OF ATTACHMENT 

We follow the same method we used to derive Eq. Taking the probability of attachment independent of in- 

and outdegree, and also identical for all nodes, 11 = l/iVx, we obtain the rate equations for the outgoing links: 

= m m lor fc > . (CI) 

Substituting UkN^ for Nk, and using dN^/dt — 1, yields equations for the density of nodes: 

iVx^^ = muk-i - (m + l)nfc . (C2) 

In the steady state, we have 

(m+ l)n^ = mn^„i . (C3) 

The solution of this equation yields the steady-state outdegree distribution: 

nl = n*c'' , (C4) 

where c = m/(l + m). Imposing the normalization condition, '"■fe = find, tIq = 1 — c = 1/(1 + m). 

We can also derive the scaling form of this distribution for comparison with the scaling form of Eq. I|ll|l in Ap- 
pendix |21 Multiplying both sides of Eq. (|C4|) by m, we get 

\ fc+i 

if^j =[(i + i/-n-^. (C5) 

In the limit m —^ oo, this equation yields the scaling function 

q{u) = mnl — exp(— w) , (C6) 

with u = k/m. q{u) is shown in Fig. [3 

As a result, the asymptotic ratio u / {k) for the outdegree distribution is equal to unity, a characteristic of the 
exponential distribution, compared to A w 0.522 in degree-dependent case. 
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involved, this order could affect the distributions, since this mechanism is degree-dependent. We use a random order of 
assembly in our simulations to avoid a possible bias: a species, i, is picked at random, and all species with a lower index 
are linked to i with a probability as described in the text. 

[26] Here, one can also use the series representation in Eq. lAlll to obtain an approximate analytical solution. However, slow 
convergence of this series makes this method very tedious, and much less accurate. 

[27] The proof of (fc) = m is very simple: (k) = J2 ki/N = K/N = Nm/N = m. 
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